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1 ABSTRACT 

In a recently suggested method for measuring the weak phase 7 in — > 
K ± (KK*)d decays, the relative strong phase Sd in D° — > K* + K~ and 
D° -> A*~A+ decays (equivalently, in D° -> A"* + A~ and D° -> K* + K~) 
plays a role. It is shown how a study of the Dalitz plot in D° — > A + A~7r° 
^cj ■ can yield information on this phase, and the size of the data sample which 

would give a useful measurement is estimated. 
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The relative strong phases for charmed particle decays obey patterns which are not 
easily anticipated from first principles but are subject to detailed experimental study, 
for example through the construction of amplitude triangles based on experimentally 
observed decay rates [UI21E1I1]- It has also been suggested jSJIHllZj that the final-state 
phase in the doubly Cabibbo-suppressed decay D° — ► K + n~ may not be the same as 
that in the Cabibbo-favored decay D° — > K~n + , even though they should be equal in 
the flavor-SU(3) limit jSj. Methods for measuring their difference have been proposed 
EH]- A Dalitz-plot method for measuring the corresponding phase difference in 
D° — > A* + 7r~ and D° — > K*~n + makes use of the interference between K* + and K*~ 
bands in D° — > As7r + 7r~ and is compatible with zero strong phase difference [TH I12j. 

Recently the question has been raised of the relative strong phase Sd between 
D° -> K* + K~ and D° -> K*- K + decays (equivalently, in D° -> A*+A" and D° -> 
K* + K~) [T^J. This phase is important in a proposed method for measuring the weak 
phase 7 in the 5 ± — ► (KK*)dK ± decays. In the present note we point out that 5d 
may be measured very directly through the interference of K* + and K*~ bands in 
D° -> A + A~7T° decays 14 . We discuss the size of present and anticipated samples 
of this final state and indicate the attainable experimental precision for 5r>. 

1 To be submitted to Phys. Rev. D. 
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We follow the notations of Ref. and define the D decay amplitudes 
A D = A(D° -> K~K* + ), A D = -> K-K* + ), 

and their ratio 



D 



r D e ldD . (2) 



The weak phase of D° — » if if* + is negligible, so the CP conjugate amplitude is 
A(D° -> K + K*~) = A D . We further define 

^ = A( J D°^K-(K + 7rV+), A^^^^if+^-TrV-)- (3) 

The amplitudes of the K* + — > if + 7r° and iT*~ — > K~tt° decays are equal. Then the 
ratio of the amplitudes in Q is 

— r D e D . (4) 



A' D A D 

Two channels of D° — > K + K~ir° go through a resonant decay of an intermediate 
if* + or if *~. They fill two bands in the Dalitz plot (see Fig. The width of these 
bands is determined by the full width V = T^»± = (50.8 ± 0.9) MeV Namely, 
the left vertical line corresponds to m 2 K+7T o = {m K *+ — T/2) 2 , while the right one 
corresponds to m 2 K+7T o = (m K *+ + T/2) 2 . Analogous expressions determine the values 
of iTL 2 K ^ n0 along the bottom and top borders of the horizontal band. For now we 
will neglect the actual Breit-Wigner distribution of event density across the bands. 
Instead, we will assume that the resonant decays are equally likely to appear near the 
central line of a band and near its borders. We will also assume that the resonant 
decays do not fall in the regions outside the two bands. We will neglect other resonant 
decays with smaller branching ratios that are not yet detected but may contribute to 
the Dalitz plot, such as D° -> tt° (K+K')^ D° -> n° (K + K~) ao , D° -> tt° (K+K~) h , 
D° — > K~ (K + 7r )^*( 14 3 )+, and D° — >■ K~ (if + 7r°) K ( 800 )+. Some of them are discussed 
later in the text and in Appendix B. Non-resonant decays uniformly fill the allowed 
phase space and provide a small background. For simplicity of the argument we will 
neglect it as well. 

The square at the intersection of the bands is the region where two channels 
interfere with each other. We denote e to be the fraction of D° — > K~ (K + tc°)k*+ 
decays that fall into the square region. This fraction only depends on masses and 
spins of particles involved in the process and the width T = T^.±. So, the probability 
of a D° K + (K~tt°) k *- decay falling into the square region is e as well. This 
probability is calculated in Appendix A: e ~ 0.039. 

Now we can write the number of decays detected in the square region of the Dalitz 
diagram: 

iV s cx \^A' D + ^A' D \ 2 = e(l + 2r D cos5 D + r 2 D )\A' D \ 2 , (5) 
while the rest of the resonant decays contribute to the bands outside the square region: 

N out oc (1 - e) (\A' D \ 2 + \A' D \ 2 ) = (1 - e) (1 + r 2 D ) \A' D \ 2 , (6) 
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Figure 1: The Dalitz plots of the D° — > K + K~tt° decay. Top panel: constructive 
interference (cosSd = 1), 113 events in the square region; bottom panel: destructive 
interference (cos 5 D = —1), 4 events in the square region. The total number of events 
in the bands is iV = 1500 in both cases. 
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so that the total number of the events detected in the bands is 



N = N S + N out oc (1 + 2er D cosS D + r 2 D ) \A! D \ 



(7) 



Experimental measurements of N s and A" provide a way of measuring the strong 
phase 5d- 

N s /N-e 



cos S 



D 



l+r 2 D 



2er D 1 - N s /N 

The uncertainty in e can be neglected because it is determined by the uncertainties in 
particles' masses and width T, which are small. The ratio td defined by Eq. (J2J) can be 
calculated from the measured branching ratios: B(D° — > K + K*~) = (2.0 ± 1.1) • 10~ 3 
and B(D° — > K~K* + ) = (3.8 ± 0.8) • 10~ 3 [T3]. Assuming the uncertainties of these 
two measurements are uncorrelated, td = 0.73 ± 0.21. These values are based on a 
sample of 35 D° —>■ KK* decays jlfij . For a larger sample, the relative uncertainty in 
td will decrease as 1/y/N. Taking the uncertainties of the decay numbers N s and A" 
to be their square roots, we can calculate the uncertainty <j(cos5d). One can show 
that the uncertainty in cos^d is mostly determined by the uncertainty in N s : 



<r(cos Sd) 



dcos5n 



a{N s 



1 



D) 



N s /N 



1 



2er D (1 - N s /N) 2 



(9) 



Unlike cos<5d itself, the uncertainty of this quantity depends not only on the ratio 
N s /N but on the total number A^ of the events detected in the bands as well. 

As an aside, note that Eq. (JHJ) predicts a linear dependence of cos<5d on Z = 
(N s /N)/(1-N s /N) with the slope S = (l-e)(l+rf,)/(2er ). We could alternatively 
write Eq. (jUJ) as 



<t(cos5d) 



d cos 5d 



dZ 



a{Z) w S 



N s /N 



1 



(1 - Ay AO 2 y/N 



(10) 



The maximum possible value of the ratio N s /N is achieved if the contributions 
from two bands are fully coherent, i.e., if cos^d = 1. In this case 



Ng_ 
N 



Ns 
N 



e(l + r D ) 2 
1 + 2er D + r 2 D 



0.074 ±0.003 



The minimum possible N s /N is a result of the fully destructive interference at cos Sd 
— 1. Then, 

N 1 = [N 1 \ _e(l-r D ) 2 

A^ 



A' 



1 - 2er D + r 2 D 



0.0020 ±0.0035 



(12) 



Thus, if cos^d is close to —1, one may observe no events in the square region. The 
source of the uncertainties in the maximum and minimum values of the N s /N ratio 
is the current 30% error in td which will be improved as more D° — > KK* decays 
are detected. Within la uncertainty, we can expect the N s /N ratio to lie between 
and 0.077. 
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Figure El shows the contours of constant a (cos 5 d) calculated for this region of 
N s /N from Eq. @ for the total number of band events N between 100 and 1500. 
The uncertainty in cosSd is an increasing function of N s /N. So, cos So will be 
measured least precisely if it is close to unity. This corresponds to a near maximum 
value of the N s /N ratio. To estimate the largest uncertainty for different numbers of 
band events, we calculate how a(cosSn) decreases with N when N s /N is fixed at its 
maximum value of 0.077: 

8.4 

CWr (cos 8 D ) PS -j= . (13) 

Now we discuss the consequences of the fact that the event density across a res- 
onant decay band is not uniform but follows the Breit-Wigner distribution. The 
differential cross-section for any point on the Dalitz plot (see Appendix A) is 



d 2 Y 

oc 



Ai(m K+7r o, m K - n o) r D e l&D A 2 (m K+n o, m K - n o) 



(14) 

The Breit-Wigner factors in the denominators make the population density nonuni- 
form across the bands while the kinematic factors Ai(mx+ n °, m K-n°) and 
A2(mK+n°, m K-n°) are responsible for a characteristic emptiness in the middle of the 
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bands. The results of a Monte Carlo simulation of the distribution ()14j) are shown in 
Fig. El 

We simulated the Dalitz plot distributions 10 times for each of 11 values of cos Sd 
between —1 and 1. For the purposes of these simulations we assumed that td is 
equal to its current central value of 0.73. The plot of Z = (N s /N)/(1 — N s /N) as 
a function of cos^d is shown in Fig. HJ To estimate a (cos 5 d) we will assume that 
the linear relationship between the two quantities still holds. Then, the slope is S = 
l/(0.0283±0.0005) = 35.3±0.6 while the maximum value of N s /N is 0.0637±0.0019 
at cos Sd = 1- Both errors are purely statistical Monte Carlo uncertainties. These 
new values of the slope S and (N s /N) max can be plugged into Eq. (fTU|) to give our 
best estimate of the maximum uncertainty in cosSd- ct(cosSd) = (10.16 ±Q.2Q)/\/N, 
with the upper bound 

0- m ax(cOsS D ) ~ — F= • (15) 

Thus, we see that the most precise measurements will be made if cos Sd is close 
to —1. The uncertainty of the least precise measurements (in case cosSd is unity) 
becomes smaller than 0.33 at N ~ 1000. Although this uncertainty is rather large, it 
at least allows one to distinguish cos^d from 0. The measurement of cos Sd will be 
improved to reach the uncertainty of 0.27 or better when 1500 resonant events are 
detected in the bands. 

In fact, 1500 resonant decays in the bands is the largest sample one can expect 
from CLEO-c. The CESR accelerator will operate at a center-of-mass energy of y/s ~ 
3.77 GeV (ip") for approximately one year. The anticipated integrated luminosity will 
reach 3 fb _1 . This corresponds to a sample of 30 million DD pairs, with 17.5 million 
of them being D°D pairs. The expected sample will exceed the Mark III experiment 
dataset by a factor of 300. Approximately 5 million of D° and D mesons will be 
flavor tagged ^7j. The other D of a pair may decay to the K + K~n° final state 
through an intermediate K*. The branching ratios of these resonant decays are 
B(D° -> K+ (K-ir°) K *~) = |(2.0 ± 1.1) • 10" 3 and B(D° -> K~ (K+ir°) K *+) = 
| (3.8 ± 0.8) • 10~ 3 , adding up to about 2 • 10~ 3 . Neglect interference effects and the 
number of decays should be around 10000. The estimated reconstruction efficiency 
for these 3-body decays is approximately 30%, so 3000 events will be detected. The 
Breit-Wigner distribution dictates that the bands of the Dalitz plot will be populated 
by half of these, i.e., by 1500 events. 

The method that will be used in data analysis will likely adopt the multi- variable 
fitting described in [TH] and ^H] instead of taking a close look at the number of 
events in the square region. We hope, however, that this note gives a good esti- 
mate of the expected uncertainty and its dependence on the total number of de- 
tected D° — > K + K*~ and D° — > K~K* + events. Other resonant decays with smaller 
branching ratios, D° -> it (K + K~)^ D° -> vr° (K + K~) ao , D° -> ir°(K + K-) fo} 
D° — > K~ (-ft' + 7T°)x*(i43o)+5 an d D° — > K~ (fC + 7r ) K ( 800 )+, may contribute to the 
Dalitz plot. The estimate of the uncertainty is most sensitive to the number of events 
inside the square region. Unless the bands of those decays overlap with it, they should 
not considerably change our estimate. 
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Figure 3: Two examples of realistic Dalitz plots of the D° — > K + K'n° decay. Top 
panel: constructive interference (cos<5d = 1), 88 events in the square region; bottom 
panel: destructive interference (cos#d = —1), 18 events in the square region. The 
total number of events in the bands is N = 1500 in both cases. 
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Figure 4: Z = (N s /N)/(1 — N s /N) as a function of cos Sd- The solid line with 
the slope of 0.0283 ± 0.0005 is the best linear fit to the results of the Monte Carlo 
simulations. The dash-dotted line is the prediction of the simplified model which 
doesn't take into account the Breit-Wigner resonant shapes (Eq. (JHJ)). 

Among the five decays listed above, only those of the k(800) + have the potential 
to contribute to the square region. However, the K is not likely to be among the 
intermediate states that make a significant contribution to D° — > K~K + tt° decays 
(see Appendix B). The meson is a narrow vector resonance which is not much 
heavier than the combined mass of two charged K mesons. Therefore, it could only 
produce a narrow diagonal band at the very edge of the Dalitz plot. Its presence 
would not change the K* band population. The same is true for a (980) and / (980) 
decays. They are lighter and broader (40 — 100 MeV) but yet not broad enough to 
significantly affect even the outer ends of the K* bands. Such a possibility is present 
for i^o(1430) decays. The square region lies outside the A"q(1430) bands and their 
impact on the number of events N s inside the square is insignificant. They can only 
make a relatively small contribution to the total number of band events N which 
would add just a small correction to the uncertainty in the strong phase 5d- 
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Figure 5: The D° — > K (K + tt°)k*+ decay in the rest frame of K + and 7r°. 
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Appendix A: Kinematics and decay amplitudes 

The first stage of the D° — > K~ (K + tt°) k *+ process is the decay of a pseudoscalar 
meson D° into a pseudoscalar K~ and a (possibly off-shell) vector K* + . After- 
wards, the latter decays into K + and tt°. From angular momentum conservation 
the helicity of K* + is 0. The corresponding polarization vector is e^*+ = e ( - A=0 - ) = 
(|Px*+|, 0, 0, EK*+)/iri K + n o. (rrix+^o is the invariant mass of K* + and the z axis is 
chosen to point in the direction of the K* + momentum Pk*+, see Fig. |5J). 

The amplitude Ai(mK+TT°, m K-ir°) of the K* + — > K + tt° decay should be Lorentz 
invariant, i.e., it should contain a product of two 4- vectors. There is only one non- 
vanishing possibility, e^.+ {j>k+ — Pit ), since the other, e K *+ (px+ + Pn a ) = £k*+Pk*+, 
is identically zero. Then the former can be written in the rest frame of K + and tt° as 

Ai(m K +n°, ™>K-n ) oc (0,0,0, 1) (E* K+ - E* ,2p* K+ ) = 2 |p^+| cosd* , (16) 

where 9* is the angle between the negative direction of the z axis and the direction 
of the K + momentum p* K+ in the rest frame of K + and 7r°. We will keep using the 
"*" subscript for quantities determined in this frame, cos 9* is given by 

cosfl* = K * ^ £ — ^ , 17 

2|P^-Iip; | 

so 

A 1 [m K +^, m K -^) oc — — , (18) 

\Pk-\ 

where 

E K~ = ( m D° - m K- -mK+7r»)/ 2m K+n° , (19) 
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£*o = {m 2 K+7T o - m 2 K+ + mlo)/2m K+7T o , (20) 
\Pk-\ = \ 1/2 {m 2 D o,m 2 K -,m 2 K+n o)/2m K + 7r o , (21) 

Ip;o| = IPk + I , (22) 
X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz . (23) 

Including the finite resonance width Tk*+ into the K* + propagator, we can write the 
amplitude of the D° — > K~ (K + tc°)k*+ decay as 

A(V —> K [K it ) K *+) oc — 2 r ' 

As for the D° — > K + (K~7T°)k*- decay, its amplitude can be derived in a similar way 
and is equal to 



A 2 {m K+n o, m K - n o] 
m 2 K _ n0 - m 2 K *- +im K ,-T 



A(D° -» K+ (K-n°) K ,-) oc r D e iS ° - 2 f^!^ _ ? (25) 



with the kinematic factor v4 2 defined as v4 2 (mx+7r0, ^K-n ) = ^i( m x-7r°) m x+7i-°)- 
The factor r^e l5D accounts for possible differences in hadronization as vector particles 
between quarks arising from the virtual W + and spectator quarks. 

Calculation of the fraction e of resonant decays that fall into 
the square region 

For the particular case of an on-shell resonant K* + we can neglect the Breit-Wigner 
denominator of Eq. (|24)1. In this case the amplitude of the D° — > K~ (K + h°)k*+ 
decay is proportional to Ai(rriK*+, mK-n )- The kinematics of the two-body D° — > 
K-K* + and K* + -> K + n° decays determine E* K . = 1.37 GeV, E* = 0.32 GeV, 
|p^_| = 1.27 GeV and |p^| = 0.29 GeV. As a result, Eq. JTZD says: 

cos9* = 1.36 (m 2 K - T0 - 1.135) , (26) 

where m 2 (K~7T°) is in GeV 2 . Thus, the amplitude of the D° — > K~ (K + tt°)k*+ 
decays is proportional to (m^-^o — 1.135). These resonant decays fill the vertical 
band in a nonuniform way: no decays happen at the middle of the band where 
m 2 K -^ — 1.135 = 0. The majority of the events will concentrate near both band ends 
where \m 2 K -^ — 1 . 135 1 is the largest. 

Now we can calculate the fraction of D° — ► K~ (K + 7t°)k*+ decays that fall into 
the square region, 

0.84 , 1.87 

e = J (x- 1.135) 2 dx / J (x - 1.135) 2 dx = 0.039 , (27) 

0.75 ' 0.40 

where (m K *+-T/2f = 0.75 GeV 2 and (m^,+ +r/2) 2 = 0.84 GeV 2 are the boundaries 
of the square region and 0.40 and 1.87 are the boundaries of the whole band. The 
latter can be derived from Eq. (|26|). 
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This simple calculation implied that the population density of the vertical band is 
constant along any cross section of the band, i.e., at a fixed iTi 2 K ^ n0 it is independent of 
variations of Tn 2 K+ ^ across the band. A more precise discussion involves a simulation 
of the interference between the Breit-Wigner resonant shapes of Eqs. f)24j) and (|25jl . 

Appendix B: Influence of scalar resonance k 

The existence of broad scalar resonances below 1 GeV has been a controversial is- 
sue for a long time [20] • A few experiments have been able to explore the pos- 
sibility of their presence as intermediate resonant states in three-body D decays. 
The modes that were studied include D + — > 7r~7r + 7r + [19 , — > 7r~7r + 7r + |21j . 
D+ -> K-TT+1T+ [22] (E791 collaboration), D° -> K^tt+tt' PU, D° -> K~ti + ti [H] 
(CLEO), and D° — > i^°i^~7r + [23j (BaBar). The first two studies obtained evidence 
for a light (478 MeV) a resonance and measured the properties of the /o(980). The 
last four might provide some information on the presence of an intermediate S-wave 
Kit resonance. Indeed, the E791 analysis of a Dalitz plot found that the best fit to the 
data is obtained allowing for the presence of an additional scalar resonance k(800)°. 
However, neither CLEO studies found evidence for k° or its isodoublet partner k + . 
The preliminary BaBar analysis saw k at the level of la which does not allow the 
confirmation of its presence. Other types of decays could also provide a glimpse of k. 
The BES collaboration found k° as an intermediate state in J/i/i -> K*(892)°K + tt~ 
decays [21], while the FOCUS collaboration studied the interference phenomena in 
D + — ► K~7T + fi + b> decays [22]. Their data can be described by K*° interference with 
either a constant amplitude or a broad spin zero resonance. 

The D° — > K^K*^ decays discussed in this note can be affected by the possible 
presence of among the intermediate states. The bands of a broad k(800) would 
cover more than 50% of the Dalitz plot, thereby interfering with the K* bands and 
affecting their population. One would expect that in this case the total branching 
ratio of D° — > K + K~7t° decays would be considerably larger than the sum of the 
D° — > KK* modes. Indeed, in D + — > K~tt + tt + an unusually high fraction (over 
90%) of decays was found to be non-resonant by previous experiments j2E]- That 
was unusual as the non-resonant (NR) contribution in three-body decays is small in 
most other cases. That was an indication of a possible broad scalar contribution 
and motivated the recent searches for it. It was found that the complex structure 
of the Dalitz plot was best explained when the k presence is assumed j22]- Then, 
intermediate decays through the kti + state account for about 50% of decays while the 
NR fraction drops to a value of 13% more characteristic of other decays. 

The present knowledge of D° — > K T K*^ decays does not reveal a similar large non- 
resonant (or broad scalar) contribution. The current data on the resonant f3E] and 
inclusive [131 12Z] decays comes from CLEO measurements. The inclusive branching 
ratio is B(D° -»• K+R-n ) = (1.24 ± 0.35) • 10~ 3 . The branching ratios of the Kn 
resonant decays are B(D° -> K + (K~7r°) K *-) = | (2.0±1.1) • 10~ 3 = (0.67±0.37)-10- 3 
and B(D° -> K~ (K + 7r ) K *+) = § (3.8±0.8)-10" 3 = (1.27±0.27)-10~ 3 . Neglecting the 
interference between these two channels (it affects just about 4% of these decays; see 



11 



Appendix A), the two branching ratios add up to (1.93 ±0.45) • 10~ 3 , consistent with 
the inclusive branching ratio within the current large uncertainties. Basically, there 
is no room for a broad scalar resonance channel. For example, it cannot negatively 
interfere with both halves of a K* band. The phase variation across it would be 
significant (~ 90°) for a K* channel and much smaller for a broad k one. If this 
channel is strong enough to cancel half the K* decays it would contribute many times 
more than that outside the K* bands. That would contradict the smallness of the 
inclusive branching ratio. Thus, we conclude that a broad scalar k, if present, could 
only comprise a small fraction of D° — > K + K~tc° decays and would not significantly 
affect the estimate of the uncertainty in the strong phase Sd- 
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